Modeling Bell's Non-resonant Cosmic Ray Instability 
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ABSTRACT 

We have studied the non-resonant streaming insta bility of charged energetic particles moving 
through a background plasma, discovered by Bell (|2004l ). We confirm his numerical results 
regarding a significant magnetic field amplification in the system. A detailed physical picture of 
the instability development and of the magnetic field evolution is given. 



Subject headings: cosmic rays- acceleration- instabilities 



1. Introduction 

T he diff usive shock acceler ation p rocess (Krym- 
sky ([19Z3); Axford et a l. (Il977f) : Bell (1978); 



Blandford and Ostriker (|l978i )) is considered as 
the principal mechanism for the production of 
the galactic cosmic rays in supernova remnants 
(SNRs). A great strength of the random magnetic 
fields that provide the scattering of energetic par- 
ticles both upstream and downstream of a super- 
nova shock is necessary for efficient acceleration. 
It was originally suggested that this may be the 
result of a gyroresonant streaming instability that 
develops due to the presence of a diffu sive stream- 
ing of accelerated particles (Bell ('1978'), Blandford 
& Ostriker (197§)). The corresponding magneto- 
hydrodynamic (MHD) waves have wavelengths of 
the order of the gyroradii of the energetic parti- 
cles. This kind of resonant cyclotron instability 
had been suggested earlier to regulate the propa - 
gation of the galactic cosmic rays (Lerche ( 1967t ). 



Wentzel (Eizi)). If the generated random fields 
are amplified up to wave amplitudes that corre- 
spond to the strength of the mean magnetic field, 
the scattering mean free path of resonant energetic 
particles will decrease to a value comparable with 
their gyroradius. This regime of diffusion is called 
Bohm diffusion and is often used to estimate the 
maximum energy of the accelerated particles. 

However, even this rather optimistic regime of 
diffusion may ensure the acceleration of cosmic ray 
protons only up to energies of about 10^ eV (La- 
gage & Cesarsky (|l983l ). Berezhko (|l996l )) if the 
magnetic field in the remnant is comparable with 
the interstellar magnetic field, which is typically 
less than 10 /zG. Significant magnetic field ampli- 
fication is necessary for acceleration up to higher 
energies. Amplification might be possible since 
quasilinear theory of the resonant streaming insta- 
bihty allows it (see e.g. McKenzie & Volk (|1982t )). 
Since this perturbation theory breaks down in the 
case of resonant interaction of particles with high 
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amplitude MHD waves, Lucek and Bell (|2000l) per- 
formed MHD simulations combined with calcula- 
tions of energetic particle trajectories and found 
that the magnetic field can be amplified to a level 
where the random field exceeds the initial mean 
field. The necessity to perform detailed trajectory 
calculations leads to strong numerical limitations 
of such simulations. 

Different qualitative treatments of magnetic 
field amplification at supernova shocks were sug- 
gested on an analytical level (sec e.g. Bell & Lucek 
(I2OOII) Ptus kin & Zirakashvih (200a), Vladimirov 
et al. (|2006f ). Amato & Blasi (20o|)). Using the 
observed synchrotron emission, field amplification 
was included on a phenomenological level in nu- 
merical solutions of the coupled gas dynamic and 
particle acceleration equations (e.g. Berezhko et 
al. ( 20 Q2), B erezhko & Volk (2004a: 2004b. ), Volk 
et al. (|2007l )). ft became clear that such ampli- 
fication can lead to particle accelaration to knee 
energies at « 3 • f 0^^ eV and ~ according to spec- 
ulative extrapolations - possibly even beyond. 

On the basis of the dispersion relation for coUi- 
sionl ess MH D waves derived by Achterberg (f 983) , 
Bell ([2OOJ) found a non-resonant streaming insta- 
bility that had been overlooked before. He argued 
that the diffusive streaming of particles acceler- 
ated at supernova shocks may be so strong that 
it modifies the dispersion relation of MHD waves 
in an essential way. Then a non-oscillatory purely 
growing MHD mode appears at scales smaller than 
the gyroradius of the particles which excite the 
instability. The growth rate of this mode can be 
larger than the growth rate of the resonant mode. 
Since the particle trajectories are only weakly de- 
flected by small-scale magnetic inhomogeneities, 
one can avoid complicated trajectory calculations. 
Only the mean c osmic-ray flux needs to be speci- 
fied. Bell ( 20041 ) performed corresponding MHD 
simulations and showed that the magnetic field 
can indeed be significantly amplified. 

In the following we investigate this non-resonant 
instability in considerably more detail. We have 
performed MHD simulations of the instability with 
very good numerical resolution and have found a 
simplified analytical description for the magnetic 
field amplification. Our results may be applied 
to different astrophysical objects, where energetic 
particles exist. 

The present work deals with the instability in 



the presence of a non-resonant energetic particle 
population streaming through the thermal plasma. 
And it simulates the nonlinear instability devel- 
opment. In a companion paper (Zirakashvih & 
Ptuskin (j2008l ). Paper II) this modeling will be 
combined with an analytical treatment of the dif- 
fusive acceleration of particles in a plane, steady 
shock with its intrinsically non-uniform energetic 
particle distribution in the precursor. 

The paper is organized as follows: The basic 
equations are given in the next section. The scat- 
tering of cosmic ray particles by small-scale ran- 
dom magnetic fields is considered in Sect. 3. The 
cosmic ray electric current is calculated in Sect. 4. 
The non-resonant instability and the MHD simula- 
tions are described in Sect. 5 and 6. The summary 
is given in the last Section. 

2. Basic equations 

We consider a system that consists of a ther- 
mal plasma and a cosmic ray gas with a negligible 
mass density. We shall treat the gas as a mag- 
netized fluid with frozen-in magnetic field B and 
induced electric field E = — [u x B]/c. Here u is 
the mass velocity (essentially equal to the thermal 
gas velocity) and c is the velocity of light. The 
electric charge density pth and electric current den- 
sity jth of the thermal plasma that determine the 
Lorentz force F = [(jth — Pthu) x B] /c may be found 
from the quasi-neutrality condition pth = ^Pcr 
and Ampere's law [V x B] = 47r(jth + jcr)/c, re- 
spectively. Here pcr and jcr are the electric charge 
density and electric current density of the cosmic 
ray gas, respectively. The Euler equation of the 
gas motion may then be written as 

fdu , \ [B X [V X B]] ^„ 
p - + (uV)uU- ^ 4. "^^ 



I 



[(jcr - PcrU) X B] , 



(1) 



where P is the thermal gas pressure. 

The evolution of the mass density p and the 
magnetic field B are governed by the continuity 
equation 

dp 



^+V(pu).0 



and Faraday's law 

dB 

'dt 



[V X [u X B]. 



(2) 



(3) 
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The equation for the plasma energy density e = 
{pu'^/2) + (B2/(87r)) + P/(7 - 1) may be written 



as: 



de „ 



7PU [B X [u X B] 

+ 



47r 



■[jcr X B]. 



(4) 



Here 7 is the adiabatic index of the gas. The 
last term in this equation describes the mechan- 
ical work produced by cosmic rays. The electric 
charge density pcr and the electric current density 
of cosmic rays jcr may be found from the momen- 
tum distribution of cosmic rays /(p, r, t). It obeys 
the equation 



5/ 
dt 



vV/ + ^[(v-u)xB]|^ = 0. (5) 
c op 



Here q is the charge of cosmic ray particles. 

The system of Eqs. (l)-(5) describes the in- 
teraction of cosmic rays and magnetized thermal 
plasma. We shall use this system in the next sec- 
tions. 

3. Scattering by the small-scale field 

The cosmic ray momentum distribution may be 
written as / = /o -I- 5f . Here /o = (/) is the 
momentum distribution averaged over the fluctu- 
ations of the magnetic field (5B and of the plasma 
velocity i5u; 5f is the fluctuation of the momen- 
tum distribution. One can use perturbation the- 
ory for the calculation of 5f when the fluctuations 
of the magnetic field are small in comparison with 
the mean field Bq. Since we are interested in the 
investigation of considerable magnetic field ampli- 
fication, the theory of perturbations with small 
magnetic field amplitude is, generally speaking, 
not applicable. Fortunately in the case of the non- 
resonant instability that is most interesting for the 
present consideration, only small-scale fields are 
generated (see also below) and the scale of the ran- 
dom field is smaller than the particle gyroradius 
in the total magnetic field. Perturbation theory 
is applicable in this case because the particles are 
only weakly defiected on the characteristic scale 
of the random field. The theory of cosmic ray dif- 
fusion in such magnetic fields was developed by 
Dolginov and Toptygin ( 1967t ). 



The equation for the fluctuation of the momen- 
tum distribution can be found from Eq. (5): 



851 
dt 



f vV<5/ + ^[(v-uo)xBo]?^ 
c op 

= -^[(v-uo)x5B]f^. 
c op 



(6) 



Here uq is the mean mass velocity and we neglect 
the velocity perturbations 5\i. It is assumed that 
they are small in comparison with the mean ve- 
locity Uq. Such perturbations would only result 
in second order Fermi acceleration which is a fac- 
tor of /c^ slower in comparison with particle 
scattering. 

Since the random field i5B has small spatial 
scales, the last term on the left-hand side of this 
equation may be also neglected. This means that 
to lowest order the particles can be considered 
to move in straight-line orbits. The calculations 
are in this case significantly simplified in compar- 
ison with the general case, when the integration 
along the helical orbits of particles results in the 
appearance of series conta ining Bessel functions 
(see e.g. Berezinskii et al. ( 199Cll )). Although the 
zero-order orbits are different for these two meth- 
ods, they give similar results in the case of almost 
any small-scale random magnetic field. The for- 
mal limit i?o ^ is not trivial and was considered 



by Tsytovich (|l97<f l. 

As long as the gradient lengths of /o and Uq 
are large compared to the spatial fiuctuation scale 
, the fiuctuation amplitudes carry a corre- 
sponding parametrical spatial dependence. We 
shall assume this in the following. Fourier trans- 
forming then in time and space we find the expres- 
sion for the Fourier transform of the distribution 
5fu^M = J d^rdt6fexp{i{ujt - kr))(27r)-4: 



<5/..k-T-^^[(v-Uo)x^B^.k]|^. (7) 
i(Lj — kv) op 



Here 5B(j ^ is the Fourier transform of the ran- 
dom magnetic field. We shall further assume that 
the magnetic field changes slowly in time in the 
frame moving with the mean mass velocity Uq and 
therefore write JB^.k (5Bk(5(aj — kup). 

Averaging now Eq. (5) and using Eq.(7) we ob- 
tain the following kinetic equation for the average 
cosmic ray distribution function /q: 
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4. Calculation of the electric current 



q Eo 



i[(v-Uo)xBo])§^ 
c J dp 



dpi dp. 



(8) 



Here Eg = —c~^ ((5u x (SB) is the mean electric 
field in the frame of reference moving with the 
mean plasma velocity Uq. 

The scattering tensor Vij in the last equation 
is determined by the spectrum of the random mag- 
netic field S,, (k) = / (i?.(r + v^)B,{v„)) exp(- 



q2 r 



X _B(r(k)(w„ - ^tom)(^;s - wos) 



(9) 



Here Ci^fc is the antisymmetric tensor. The scat- 
tering tensor makes the cosmic ray distribution 
isotropic in the frame moving with the velocity 
Up. Expression (9) may be simplified in the case 
of an isotropic random magnetic field i?ij(k) = 
Bisot r{k)(5i, - kikj/P) (Dolginov & Toptygin 



Il967l) 1 



: P^vip) 



(v - Uo)'^Sij - {vj ~ uoi){vj ~ upj) 
2v\v - uo| 



(10) 

where the scattering frequency i>{p) is given by the 
formula 



v{p) 



4 p^c^ 



d-^kBu 



.{k)/k 



(11) 



Here the spectrum of the isotropic magnetic field 
is normalized as {SB^') — J (PkBiBotr{k)■ 
ln the diffusion approximation the average cos- 
mic ray distribution function /o may be written 
as 



/o(p)=iVb) 



— pJ(p), 

pv 



(12) 



where N{p) is the isotropic part and J(p) is the 
cosmic ray flux density, which is the sum of the 
diffusive and advective flux densities and is given 
by Eq. (A.3) in Appendix A. 



We shall use Eq. (7) for the calculation of the 
electric current density of the cosmic ray gas. Sub- 
stitution of the expression (12) into Eq. (7), mul- 
tiplication by gv, and integration over momentum 
space give the fluctuating part of the Fourier trans- 
form of the cosmic ray electric current (5jcr,cj,k- 
Disregarding terms of the order uq jv we obtain 



ikr): 



37r^5(a)-kuo) / — <5(kv)(J Jv x ,5B 

c J pv 

(13) 



k V. 



Here Jd — J -I- Uq | ^ is the average cosmic ray 
diffusion flux. The appearance of the (5-function 
in this equation is due to the Landau resonance 
a; = kv (cf Lifshitz and Pitaevskii 1 198l[ )') in 
Eq.(7). 

The total electric current of the cosmic ray gas 
is the flux J multiplied by the particle charge q. 

Performing the integration on the two angles 
in momentum space, and calculating the inverse 
Fourier transformation of the last equation, we ob- 
tain the expression for the diffusive electric current 
jd — jcr — PcrU that appears in the right-hand side 
ofEq. (1): 

jd = g y" d^P^d cxp(ik(r - M^t)) 

X j (fc'[<5Bk X J,] - k(k[5Bk X J,])) . (14) 

The flrst term on the right-hand-side of this 
equation is simply the zero-order term of the ex- 
pansion in the magnetic fluctuation i5B, while the 
second term is the linear term of the expansion. 
The latter is always smaller than the first term if 
the small-scale field approximation is valid, that 
is if qB/pck << 1, which means that the particle 
gyroradius pc/qB is large compared to the scale 
k^^ of the magnetic field. 

However, in some cases the second integral may 
play a role. When the cosmic ray streaming is not 
strong, it results only in a small change of the 
dispersion relation of MHD waves. The first in- 
tegral in Eq. (14) in this case produces only a 
small shift of the frequency of MHD waves. The 
second term then gives a small imaginary part of 
the frequency and describes a resonant wave insta- 
bility based on the Landau resonance in Eq. (7). 
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Within the hmits of the small-scale field approxi- 
mation, and for strong cosmic ray streaming, this 
resonant instability is ineffective in comparison to 
the well-known gyroresonant streaming instability. 

Another important point is that the second in- 
tegral on the right-hand side of Eq. (14) ap- 
pears in the calculation of the me an for ce act- 
ing on the thermal plasma (Ptuskin ( 1984( l): F = 

C"M[B X (j„. -pcrUo)]). 

Using Eq. (14) and averaging we obtain 



-[jdXB]c-i 



F = 

4 c2 



^Jd / d^kBisotrik)/k 



(15) 



It was assumed here that the random field is 
isotropic and expression (11) was used. Since the 



diffusive fiux is equal to 3d 
fusion approximation, F = 



-VP,. 



VN in the dif- 

r, where Per = 



J d'^ppvN(p)/3 is the cosmic ray pressure. In the 
sense of MHD theory, Eq. (1) then describes the 
overall momentum balance of the system, where 
the forces on the r.h.s. are the Lorentz force and 
the gradient of the overall pressure, thermal plus 
nonthermal. Indeed c"^ ([(jcr ~ PcrUo) x B]) = 
— VPcr is the general form of the average cosmic 
ray momentum balance. 

In the sequel we will n eglect the 2nd term of Eq. 
(14) (see also Bell (|2004l) ). This is justified in the 
linear analysis and in the nonlinear simulation of 
the non-resonant instability if the small-scale field 
approximation is valid. 

5. Non-resonant streaming instability 

The dispersion relation for small-scale MHD 
perturbations may be found from Eqs. (l)-(5). 
In the case when the cosmic ray diffusion flux Jd 
and the wavenumber k are parallel to the mean 
magnetic field B n this dispersion relation may be 
written as (Bell |2004h ): 



jdBo 
/ 

cpo 



(16) 



Here po is the mean plasma density, — 
Bo I V^Ti'Po is the Alfven velocity, jd = 9 / rf'^pJd 
is the average diffusive electric current of cos- 
mic rays, and the two signs ^ correspond to 
the two circular polarizations. A non-resonantly 




Fig. 1. — Explanation of the non-resonant stream- 
ing instability. The magnetic spiral (thin solid 
line) is stretched by the Lorentz force F = — [jd x 
B]c~^ that appears due to the diffusive cosmic ray 
electric current jd. 



unstable MHD mode appears if the condition 
k < kc — jiiBo/cpoV^ is fulfilled. The unsta- 
ble magnetic field line spiral expands in the di- 
rection perpendicular to the mean magnetic field 
(see Fig.l). The mode with k = fcc/2 has the 
maximum growth rate 7niax ■ 

_ idBp 

7max — „ , , J ) 
2cpoVa 

which does not depend on the magnetic field 
strength. 

Since we assumed that the scale oi the per- 
turbations is smaller than the gyroradius pc/qBo 
of the energetic particles, the wavenumber k 
should obey the condition qBo/pc << kc- This 
means that the necessary condition for instability 
is kc > qBo/pc. This condition may be rewritten 
as 



>> 1 



(18) 



V B^/Att 

Here Wcr arid Ccr are the bulk velocity and the en- 
ergy density of the cosmic ray gas, respectively, 
and V is the velocity of energetic particles. This 
condition is easily fulfilled at the shocks of SNRs, 
where u^r is of the order of the shock velocity ui 
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and the energy density of the relativistic particles 
may be comparable with pw^. 



6. Numerical modeling of the non-resonant 
instability 

We have numerically modeled the non-resonant 
instability similar in spirit to the modeling of Bell 
(|2004l ). The MHD Eqs. (l)-(4), written in dimen- 
sionless form, were solved numeric ally. W e used 
the numerical method of Pen et al. (j2003[ ) . It is a 
second order in space and time, flux-conservative 
total variation diminishing MHD scheme which 
enforces the VB — constraint to machine pre- 
cision. The nonlinear flux limitcr "minmod" was 
used. 

The dimensionless time t, the space coordinate 
z and the velocity u are defined as t = tVako, 
z = koz, u = u/Va, respectively. Here fco is the 
wavenumber that corresponds to the real size of 
the numerical box 2n/kQ. The dimensionless den- 
sity p and the electric current J can be expressed 
via the magnetic field i?o and the Alfven velocity 
14 as p = ^TTpV^/Bl and J = iirj/ckoBo. The 
dimensionless wavenumber k is simply k = k/k^. 

The simulations were performed in a cubic box 
with size 27r. Periodic boundary conditions were 
imposed at the sides of the box. We used 256'^ 
grid cells in our simulations (to be compared with 
128^ cells in BeU's case). 

At t = the plasma pressure and density are 
uniformly distributed in space. Small random 
magnetic perturbations corresponding to isotrop- 
ically distributed Alfven waves with a one dimen- 



1/2 

sional spectrum (x k~^ and (-B^) / Bq 
were added to the mean unit strength magnetic 
field that is in z direction. Here () denote the spa- 
tial average over the simulation volume. We use 
the values 7 = 5/3 and /3 = 1. Here (3 = AirP/B^. 

The evolution of the magnetic field and the 
mass velocity fluctuations, together with the evo- 
lution of the sound speed, the mean electric field 
Eq and the characteristic scale Lm of the mag- 
netic field, for a dimensionless cosmic ray current 
J = 16, is shown in Fig. 2. The characteristic scale 
Lm = J dk,B±{k,)\k,\-'^/ J dk.BjLik,) is deter- 
mined via the spectrum B±{kz) of the perpen- 
dicular component of the random magnetic field 
B±ikz) = J dk^dky(B^^Ck) + Byy{k)). 

After a brief initial stage the fiuctuations grow 



0.09 



100 



-(£o/So)^/'(c/K)i/> 




Fig. 2. — Numerical results of the modeling of 
the non-resonant instability with the dimension- 
less cosmic ray electric current J = 16 as a func- 
tion of normalized time t. The r.m.s. values of 

1/2 

the magnetic field fluctuation {SB'^') divided 
by the mean magnetic field strength Bq and of 

1 /2 

the mass velocity fluctuations (^u^) divided by 
the Alfven velocity Va are shown by the solid and 
dashed lines, respectively. The ratio of the sonic 
velocity Cs and the Alfven velocity Va is shown by 
the dotted line. The square root of the ratio of 
the mean electric field Eq and the mean magnetic 
field Bq multiplied on (c/Va)^^^ is shown by the 
thin solid line. The characteristic scale i,„ of the 
magnetic field is also shown by the thin dashed 
line. 
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t % 

Fig. 3. — Slices of the magnetic field strength 
through the center of the box, obtained at i = 0.82 
(top), I = 1.03 (middle) and I = 1.48 (bottom), in 
the run with J = 16. The scaling of the magnetic 

1 /2 

field strength is logarithmic between 0.32 (-B^) 

1 /2 

(white) and 3.2(5^) (black). The mean mag- 
netic field and the diffusive electric current are in 
the z direction. 




Fig. 4. — The one-dimcnsional spectra of the per- 
pendicular component of the random magnetic 
field B±{kz) obtained at f = (thin solid line), 
i = 0.54 (sohd line), i = 0.74 (dashed line), 
i = 0.97 (dotted line) and i = 1.26 (thick solid 
line). All spectra are normalized to the square of 
the mean magnetic field. The critical wavenumber 
fee = 16 is marked on the x-axis. 

exponentially with a growth rate that is slightly 
smaller than 7max = 8 in dimensionless units. The 
initial growth of the magnetic fluctuations is not 
exponential, since only a part of the initial pertur- 
bation corresponds to unstable modes. 

At t = 0.7 the magnetic perturbations are al- 
ready comparable with the mean magnetic field. 
Parts of the magnetic spiral expanding into the 
Xy-planc begin to collide with their surroimdings. 
When this happens multiple shocks are formed. 
The shape of these shocks may be seen in the top 
panel of Fig. 3, where the magnetic field strength 
in perpendicular slices through the center of the 
box are shown. In the YZ-plane, which contains 
the initial magnetic field vector, they look like bow 
shocks. The shocks corresponding to the adjacent 
turns of the same magnetic spiral are clearly seen. 
These shocks are almost circular in the XY-planc. 
Low density cavities appear inside these shocks. 
The size of these cavities in the XY plane is larger 
than the size in the z direction. At later times the 
shocks collide with each other in the XY plane 
and the gas motion becomes strongly turbulent 
(middle and bottom panels of Fig.3). 

The MHD turbulence has rather small spatial 
scales in the early stage of the magnetic field 
growth. The scale of the magnetic field increases 
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with time (see Fig. 2). This increase is slower 
in comparison with what was suggested by Bell 
k ^ Anjd/Bc (or k ^ JBq/B in dimensionless 
units) when the magnetic tension forces are com- 
parable with the Lorentz force produced by the 
cosmic ray electric current. This effect is illus- 
trated in Fig. 4, where the magnetic spectra ob- 
tained for several instants of time, are shown. It 
is clear from this figure that a nonlinear trans- 
fer of magnetic energy takes place in the system. 
The increase of the magnetic energy at wavenum- 
bers k > kc that are not excited in the linear ap- 
proximation demonstrates the non-linear transfer 
of energy to smaller scales. The magnetic energy 
in the small wavenumbers also grows faster than 
predicted by the analytical growth rate formula. 
For example the amplitude of the harmonics with 
wavenumber k — 1 increases by a factor of 4.4 
during the period from i = 0.74 up to f = 0.97 
(see Fig. 4). The linear growth rate for this har- 
monics 7 = \/ Jk — k"^ = \/T5 corresponds to an 
amplification factor of about 2.2 during this pe- 
riod. This demonstrates the nonlinear transfer of 
energy to larger scales. Thus the magnetic en- 
ergy is non-linearly transferred to both smaller 
and larger scales compared to the linearly excited 
spatial scales. 

Towards the end of the simulation, at t ~ 1.5, 
the scale of the magnetic field is comparable with 
the size of the box. At this point in time the in- 
ternal energy of the thermal gas is roughly equal 
to the kinetic and magnetic energy in our simu- 
lation. Bell continued his calculation beyond this 
point and found that at later times the magnetic 
field reaches a saturation value. Also continuing 
our simulations we found that the magnetic field 
continues to grow. We believe that this difference 
is due to the fact that different MHD codes were 
used. The difference is not important however, 
since the simulation does not model the real sit- 
uation any more when the scale of the field has 
become comparable with the size of the spatial 
simulation domain. In this sense our results are 
similar to Bell's results. 

The growth of the MHD perturbations de- 
creases with increasing magnetic field amplifi- 
cation. As estimated in Appendix B, the mag- 
netic field is amplified only linearly at large times, 
B / Bq ~ Jt in dimensionless units, whereas the 
gas thermal energy density increases oc t^. These 



dependencies are derived using the equation for 
the evolution of the magnetic helicity. 

7. Conclusion 

We have modeled the non-resonant instability 
produced by a flux of charged energetic parti- 
cles that is driven through a scattering thermal 
plasma. Using a significantly better numerical res- 
olution we basi cally confirm the results obtained 
earlier by Bell ( 20041 ). The magnetic field may be 
amplified significantly. The unstable magnetic spi- 
rals collide with each other and shocks of moder- 
ate strength are formed as the instability develops 
(see previous Section). These shocks lead to sig- 
nificant gas heating. Since free expansion of the 
magnetic spirals after collision is impossible, the 
field grows only linearly in time at later epochs 
(see Eq. (B3)). If the system has enough time to 
evolve, the magnetic field growth will be stopped 
when the gyroradius pc/qB of the energetic par- 
ticles in the amplified field B will drop down to 
the scale of the amplified field k^^ ^ cB/Anj^. 
This determines the v alue of the saturated mag- 
netic field (Bell (|2004f ). Pelletier et al. (|2006l )): 



B^ 
An 



(19) 



We should note that the small-scale approxima- 
tion considered in Sect. 3 becomes invalid when 
the magnetic field reaches this saturation value. 

Since the instability is driven by the Lorentz 
force, the corresponding MHD turbulence has spe- 
cific properties. It has non-zero magnetic helicity 
and a non-zero mean electric field Eq parallel to 
the mean magnetic field. 

We expect that the formation of multiple shocks 
in three-dimensional MHD turbulence will also oc- 
cur for other instabilities, in particular for the res- 
onant streaming instability, driven by cosmic rays. 

The scattering of energetic particles by the 
small-scale magnetic inhomogeneities can be de- 
scribed using the Dolginov- Topty gin approxima- 
tion (Dolginov & Toptygin (jl967t )). The appear- 
ance of a mean second order electric field Eq = 
— c^^ ((5u X (5B), which is oppositely directed to 
the electric current of the energetic particles, mod- 
ifies the cosmic ray transport equation (see Ap- 
pendix A). 

This non-resonant streaming instability may 
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be important in any astrophysical site where a 
strong electric current of energetic particles exists 
and where the initial magnetic strength is small 
enough (see condition (18)). Supernova remnants, 
star burst galaxies, galaxy cluster accretion shocks 
and AGN jets are possible candidates for an ap- 
plication of this instability. 

The results obtained will be used in_a compan- 
ion paper by Zirakashvili & Ptuskin ( 20081 ) (Paper 
II) for a model of diffusive shock acceleration in 
young SNRs in the presence of the non-resonant 
streaming instability. 



We thank the anonymous referee for a number 
of valuable comments. VSP and VNZ acknowl- 
edge the hospitality of the Max-Planck-Institut fiir 
Kernphysik, where this work was mainly carried 
out. The work was also supported by the RFBR 
grant in Troitsk. 
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A. Diffusion approximation in the presence of the additional electric field 

The cosmic ray transport equat ion is modified in the presence of the additional mean electric field Eg = 
-c-i (<5u X SB) (Fedorov et al. (|l992h '). Let us substitute the cosmic ray momentum distribution (12) 
into Eq. (8) with the scattering tensor (10). Performing the expansion up to the second order in u/v and 
collecting the terms independent of the direction of the particle velocity v and separately those proportional 
to the velocity we obtain after some algebra 



ON 
'dt 

dJ 



VJ 



1 d p2 



-q(EoJ) = ^ — J 

op V op 



U07 



pdN 
3 9^ 



yVA^ + .Eo-- 



J + Uq 



pdN 
3 dp 



fl X J + Uq 



[fl X J] 

pdN 



3 dp 



(Al) 



(A2) 



Here = guBo/pc. Assuming a slow time evolution we neglect the time derivative in the last equation. 
Then it can be used to find the cosmic ray flux J: 



Oj 



qdN 
V dp 



pdN 
3 dp ' 



(A3) 



where the diffusion tensor Dij has the following form 



= (^11 - D±_)bihj + D^6.,j + DAEijkbk- 



(A4) 



Here b = B0/-B0 is the unit vector in the direction of the mean field Bq, -D|| and D± denote the parallel and 
perpendicular diffusion coefficie nts, re spectively, and Da is the antisymmetric diffusion coefficient. They are 
given by Dolginov & Toptygin ( 1967t )): 



V 

3^' 



2 " 



(A5) 



Then Eq. (Al) reduces to 



dN 



pdN 



— + uoViV - ^^Vuo = V, + — — i^E, 



dt 



3 dp 



1 d p^q 



p2 Qp u 



D, 



V dp 



(A6) 



Eq. (A6) shows how the presence of the mean electric field Eq modifies the diffusion term of the cosmic ray 
transport equation. 



B. Magnetic helicity 

As noted first by Pelletier et al. ( 20061 ). the MHD turbulence generated by the non-resonant instability 



has nonzero magnetic helicity H = (SASB) , where SA is the perturbation of the magnetic potential. The 
magnetic helicity iJ is a us eful qu antity and it is often used in the theory of MHD turbulence and dynamo 
theories (see e.g. Biskamp ( 20031 )). Faraday's equation (3) may be used for the determinatio n of th e time 
evolution of this quantity. For the periodic system we have (cf. Subramanian & Brandenburg (|2004l) 1 



dH 
'dt 



^2cBoEo 



(Bl) 



We neglect magnetic dissipation here. This seems well justified because the magnetic helicity is an integral 
quantity and is not transferred by nonlinear interactions to smaller scales where dissipation is essential. In 
this sense the magnetic helicity is different from the nonthermal (kinetic + magnetic) energy that may be 
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transferred to smaller and smaller scales where it is transformed into gas thermal energy even in the case of 
infinitely small viscosity 

The mean electric field appears in the system as a response of the medium to the cosmic ray electric 
current. The evolution of the energy density of the plasma takes place according to Eq. (4). The right-hand 
side of this equation is simply — Eojd- Comparing with Eq. (Bl) for a time independent jd, we obtain the 
relation 

H - = const (t) (B2) 

Jd 

It is worth emphasizing that Eq. (B2) is valid for the total plasma energy density e which appears in Eq. (4). 
During the stage of exponential growth the instability produces magnetic field and velocity perturbations, 

whereas gas heating is important at later times. 

We may use Eq. (Bl) for a derivation of the equation for the magnetic field amplification. Since the 
turbulence is helical, the magnetic helicity (that is the product of the magnetic field and the vector magnetic 
potential) is _ff ^ B'^/k, the electric field is Eq ^ V^B'^ /cBq, and k Anjd/Bc. Then Eq. (B.l) yields the 
equation for the amplification of the magnetic field: 

m- ^ ^^^^ 

The numerical factor in this equation is of order unity, according to our numerical results. 

There is another way to obtain this last equation. The field is amplified by turbulent motions of the 

medium, that is dB/dt ^ B^kuk, where u^. is the turbulent velocity with wavenumber k. Assuming equipar- 
tition, Uk ~ Vq^B/Bq, and the estimate kB ~ 47rjd/c, we arrive at Eq. (B3). 

Since H oc B^ at late times, also (e) oc t^, cf. Eq. (B2). On the other hand B'^ oc t^, and therefore at 
very late times the gas internal energy dominates. 
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